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Abstract 

Iron abundance distribution is now known for 12 clusters of galax- 
ies. For some clusters (e.g. Centaurus) the observed abundance in- 
creases toward the cluster center, while for the others (e.g. Coma and 
Hydra-A) no significant inhomogeneity was observed. In order to un- 
derstand this difference, we investigate the influence of cooling flow 
and turbulence produced by galactic motion on the iron abundance 
distribution by simple spherical models. We show that the cooling flow 
has a significant effect to flatten the iron abundance distribution if the 
flow velocity is sufficiently large. Further, by applying our analysis to 
the above clusters we show that we can give a systematic account for 
the observed variety of the iron abundance distribution qualitatively. 



1 INTRODUCTION 



It is now well established by X-ray observations that intracluster 
gas is contaminated by heavy elements up to the level of 1/4 to 1/2 
of the solar abundance. However, until recently we have had little 
information on the spatial distribution of heavy elements inside each 
cluster of galaxies due to the poor spatial resolution of detectors. The 
exceptional cases were the Virgo, the Coma and the Perseus clus- 
ters. For the Virgo it was found that the iron abundance increases 
toward its center (Koyama, Takano, & Tawara 1991). In contrast for 
the Coma no significant inhomogeneity of the iron abundance was ob- 
served (Hughes, Gorenstein, & Fabricant 1988; Hughes et al. 1993). 
For the Perseus inhomogeneity of the iron abundance is a matter of 
argument (Ponman et al. 1990; Kowalski et al. 1993; Ohashi et al. 
1994). 

Some explanations for this difference in the iron abundance dis- 
tribution have been presented so far. For the Virgo cluster Okazaki 
and others showed that the central concentration of iron is explained 
if one assumes that iron was ejected only from early type galaxies 
and that the iron distribution has not changed after the initial star 
burst (Okazaki et al. 1993). On the other hand for the Coma cluster 
Fabian et al. (1994) argued that its double-center structure suggests 
some violent merging phenomenon in its past, which smeared out the 
inhomogeneity of the iron distribution produced at its early phase. 

The recent observation by ASCA (Ikebe et al. 1994 ; Fukazawa 
1994 ; Fukazawa et al. 1994 ; Ohashi 1994) and the combined analysis 
of the Ginga LAC and Einstein SSS (White et al. 1994) have changed 
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the situation of this problem significantly. We now have data on the 
spatial distribution of iron for 9 more clusters; Centaurus, AWM7, 
A496, A2142, A2199, A1060, A1795, MKW3s, and Hydra-A. The first 
five of them as well as Virgo have centrally concentrated distributions 
of iron, while for the rest no significant inhomogeneity was observed. 
(However White et al. show that the possibility of a weak iron con- 
centration in A1795 cannot be completely discarded.) Though A1060 
has the double-centered structure like Coma, the other clusters of the 
latter group has no such structure. So the above explanation for Coma 
does not apply to the latter group. Instead they have strong cooling 
flows. This suggests some relation between cooling flow and the iron 
abundance distribution. 

Motivated by these observations, in the present paper, we investi- 
gate the influence of cooling flow on the iron abundance distribution 
by simple spherical models. We also examine the influence of the tur- 
bulent mixing provoked by the motion of galaxies by assuming that it 
can be effectively treated as diffusion. 

The paper is organized as follows. First in the next section we 
describe the fundamental assumptions in detail and write down the 
basic equation to determine the time evolution of the iron abundance 
distribution. Then in §3 we solve this equation numerically for some 
typical models and quantitatively estimate the effect of cooling flow 
and motion of galaxies on the iron abundance distribution. On the 
basis of this general analysis, in §4, we construct specific models for 
each of the above 12 clusters taking into account the physical prop- 
erties of the individual clusters, and by comparing the model predic- 
tions with observational data we discuss the possibility of accounting 
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for the observed variety of the iron abundance distribution system- 
atically. Section 5 is devoted to conclusion. Throughout this paper 
Ho = 50km _1 s _1 Mpc _1 is assumed. 

2 Assumptions and Basic Equations 

We assume that the structure of a cluster of galaxies is described 
by spherically symmetric models. As a cooling flow affects only the 
central region of two or three core radii, we think that this approxima- 
tion is good. Since the atomic diffusion of iron in the intracluster gas 
is negligible in the central region, under this assumption of spherical 
symmetry the time evolution of the iron abundance Z(r, t) is described 
by the differential equation 



velocity, D is an effective diffusion coefficient to be explained later, 
and S(r,t) is the source term. 

We consider two models as the source term in this paper. In the 
first model it is assumed that iron was first produced by star bursts in 
each galaxy at its formation and ejected into the intracluster space by 
galactic winds during some initial short period < / < ^gw(< IGyrs) 
approximately at a constant rate. Further it is assumed that iron was 
ejected only from early type galaxies in proportion to their masses 
(see Arnaud et al. 1992). Under these assumptions the source term is 
expressed as 




(1) 



where r is the distance from the cluster center, v is the inward radial 
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during the windy period and vanishes after that, where p ga ± is the 
mass density of galaxies averaged over the dynamical time scale of the 
cluster, p gas is the intracluster gas density, m p is the proton mass, 
TOFe is the atomic mass of iron, /e^) is the local fraction of early-type 
galaxies in a cluster, and f w and fp e are the luminosity averaged mass 
fractions of gas and iron ejected from galaxies over the Hubble time, 
respectively. 

In the cluster centers, however, it is very likely that additional iron 
was supplied from the big central cD galaxies and galaxies passing 
nearby at their formation. In the second model this effect is taken 
into account by modifying the above source model around the center 
in the following way. 

First for simplicity we assume that the central cD galaxy affected 
the iron distribution only during the windy period, although it is at 
present controversial whether cD galaxies were formed within a short 
time after the formation of clusters and whether they affected other 
galaxies only during an early epoch (Malthmuth & Richstone 1984; 
Merritt 1984, 1985; Tonry 1987). Then, from the mass conservation, 
the metal rich gas ejected from the central cD galaxy and galaxies 
disrupted by it expands after the windy period to the region of the 
radius R cF) determined by 



where p cF) is the present stellar distribution of the cD galaxy. The 
iron abundance in this region is given by 




(3) 



/Fe m p 

fw m Fe ' 



(4) 
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Hence in the lowest approximation the source term for clusters with 
cD galaxies at center is given by 

S(r,t) = -*—Z cD (5) 
tow 

for < r < i?cD 5 an( i by Eq.(2) outside. This is the second source 
model. 

The spatial behavior of S in the first model is determined by the 
mass distribution functions p ga i, /e an( i Pgas- In the general analysis 
in §3 we assume that p ga ^ is given by the King model 

/> S al(r) = /> ga i(0)(l + r 2 /a 2 r 3/2 , (6) 
and p gas by the /3-model 

/Wo = /wo)(i + r 2 /« 2 r 3/3/2 , (7) 

and treat p ga i(0), p ga s(0) and the common core radius a as free model 
parameters. As the fraction of elliptical galaxies two cases of /e = 1 
and /e = (e~ r l a + l)/2 are considered. 

On the other hand, when we construct specific models of individual 
clusters in §4, we determine these parameters from observational data 
as far as possible. Further for some clusters different forms are adopted 
for p gas from the observational ground. The details will be given in 
that section. 

Throughout the present paper we assume that these density dis- 
tributions do not change after the windy period. This assumption will 
be good for a region within two or three core radii from the center, 
provided that there occured no violent phenomena such as merging. 

In the second model we must further specify p c ^. Observations 
show that the surface brightness profiles of the halos of cD galaxies 

6 



can be described by r~ 2 (Oemler 1976). This is consistent with the 
result of the theoretical study by Duncan, Farouki and Shapiro (1983), 
who showed that mergers rapidly establish the r~ 2 surface brightness 
profile outside the core of the cD galaxy. So it is natural to assume 
that the present star distribution of the cD galaxy, PcD(r), is given by 

/ , cD (r)= / , cD (0)(l + r 2 /a c 2 D )- 3/2 . (8) 

The inhomogeneity in the iron distribution produced by the source 
term is modified by the cooling flow, which is represented by the first 
term on the right-hand side of equation (1). As is clear from this 
equation, they affect Z only through the velocity of the cooling flow 
v. 

In the cooling flow model v(r) is related to the inward mass flow 
rate through the sphere of radius r, M(r), by 

M(r) = 4irr 2 p gas (r)v(r). (9) 

Using a multi-phase model, Thomas, Fabian, & Nulsen (1987) deter- 
mined M(r) from observations and found that M(r) increases approx- 
imately in proportion to r until around the cooling radius r c . Here 
the cooling radius is defined by the condition that the cooling time of 
the gas is equal to the Hubble time at r = r c . They also found that 
M vanishes at some finite radius from the center for some clusters. 
Though this feature is not definite enough due to the poor spatial res- 
olution of observations, theoretical arguments suggest that it is very 
likely to be real (Fabian & Nulsen 1977; Cowie, Fabian, & Nulsen 
1980, Loewenstein & Fabian 1990). Hence, in the present paper, we 
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assume that M is expressed in the form 

f for < r < r st , 

M(r) = I M c (r - r st )/(r c - r st ) for r st < r < r c , (10) 
{ M c for r c < r , 

where M c = M(r c ) and r st are the parameter. We call r st the stagna- 
tion radius. 

Since we assume that Pgas(r) is time-dependent, the decrease of 
M(r) toward the center means that a substantial amount of hot gas 
is steadily lost from the cooling flow, which is balanced with a gas 
inflow from outer region. The most probable mechanism for this is 
the thermal instability in the flows (Nulsen 1986; Malagoli, Rosner, 
k Bongo 1987; White k Sarazin 1987; Balbus 1988; Balbus k Soker 
1989). Since the thermal instability does not transfer heavy elements 
between hot regions and cool regions, we assume that the hot gas is 
removed from the cooling flow retaining its abundance. 

We further assume that the cooling flow is time-independent and 
present after winds had stopped (t > ^gw)- Since the time scale of the 
winds is < 10 9 yr (David, Forman, k Jones 1990, 1991a) and that of 
the cooling flow is ~ 0.2 - 2 X 10 10 yr (Edge, Stewart k Fabian 1992), 
this assumption seems to be reasonable. 

In contrast to the bulk inflow motion, turbulence of hot gas pro- 
duced by motion of galaxies and the associated diffusion or mixing 
of heavy elements are very intricate phenomena, and little theoretical 
investigation has been done on them. So we replace them by the fol- 
lowing simple diffusion model in this paper. Each galaxy is assumed 
to disturb a spherical region of radius R e ff, which we call the effective 
radius of the galaxy, and homogenizes the iron abundance inside it. 
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For any given point the abundance information diffuse over distance 
of order i? e ff? every time some galaxy passes by it within the distance 
i? e ff from it. Hence, assuming that the motions of galaxies are random 
and isotropic locally, the effective diffusion coefficient D is given by 

D ~ R 2 eS X n^airRlx/G = iran^R^/Q, (11) 

where a and ra ga i(r) are the local virial velocity and the local number 
density of galaxies. We assume that the r dependence of ra ga i(r) is 
given by the King model Eq.(6) with p ga i(0) replaced by n ga [(0). If 
the galactic winds are not so strong, R e ^ is expected to be at most of 
the order of the tidal radius of the cluster gravitational field, 

R eS < 60(a/500kpc)kpc (12) 

(Malthmuth k Richstone 1984; Merritt 1984). 

3 GENERAL NUMERICAL ANAL- 
YSIS 

The generic model explained in the previous section has 15 param- 
eters apart from the freedom in the choice of /e; ^gw 5 /Fe, /w> Pgal(O), 
Pgas(O), a, /3, Pcd(O), a c D, M c , r st , r c , a, n ga i(0) and R eS . In order to 
determine the present value of Z by integrating the evolution equa- 
tion, we must also specify the time from the beginning of the windy 
period, to. However, these parameters do not affect the final value of 
Z independently. In fact, if we introduce the dimensionless quantities 

T = t/t , (13) 

x = r/a, (14) 
9 



v = vto/a, (15) 
D = Dt /a 2 , (16) 
S = St , (17) 

the evolution equation (1) is written only in terms of them as 

dz „dz i d ( ~ 2 dz\ ~ 

OT OX X z OX V OX J 

Here the a;- dependence of the dimensionless quantities v, D and S are 
given by 

[ v a x- 2 (l + x 2 fM 2 7 rst/a 7 (r st /a<x <r c /a)„ nS 



v„x 



[1 + x 2 ) 3 ^ 2 (r c /a<x) 



b(x) = b (l + x 2 )- 3 ^ 2 , (20) 

S(x ' T) = T^ x { Z ztux)(i + x 2 f^)i 2 (Id/^) ,(21) 



where 



2.8 Xl0 -2 



1.33 X 10 10 yr lOOM /yr 



D = 7.3 X 10 -4 



10 2 m p /cm 3 / \250kpc, 

to o a 3 ra ga i(0) 



(22) 



1.33 X 10 10 yr 10 3 km/s 10 



V250kpc7 VlOkpc; ' 



(23) 



Zo = Z cD ^-^\. (24) 

From these equations we see that only 11 of the 16 parameters are 
independent in effect. 

In this section we reduce the number of these free parameters by 
giving typical fixed values to some of them, because the main purpose 
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of this section is to look into the characteristic features of the influence 
of the cooling flow and the galactic motion on the iron distribution 
evolution and get rough estimates of it. 

First we take a = r c = 200kpc, to = 10 10 yr, /gw = 10 8 yr, and 
p gas = 5 X 10 _3 TOp/cm 3 . For this choice the dimensionless parameter 
v a is controlled by M c . 

Next in order to fix the source term we take f w = 0.18 and 
fp e = 6.8 X 10~ 4 , which correspond to the IMF exponent x = 1.5 
(David et al. 1991b). Although David et al. (1991b) showed from the 
comparison with observations that x = 1.0 is more preferable than 
x = 1.5, Okazaki et al. (1993) pointed out that the ejected iron masses 
in Davids' models are rather overestimated because the iron which is 
used in star formation in stars is neglected in their models. Therefore 
we adopt the exponent x = 1.5. For this choice Z c d = 6.8 X 10~ 5 , 
which is 1.7 times the solar iron abundance 4 X 10~ 5 . To determine 
/°gal(0) we use the empirical relation between the number of the galax- 
ies within 500kpc from the cluster center, iVsoo, and the total luminos- 
ity of them. For M/L = 8Mq/Lq this relation is expressed in terms 
of the total mass of galaxies in the same region, M500, and iVsoo as 

M 500 = 1-44 x IO^N^Mq. (25) 

(Edge k Stewart 1991b). With this relation and Eq. (6) p ga i(0) 
is determined. We tentatively take iVsoo = 21 (Bahcall 1981) or 
Pgal(0)/rap = 6.1 X 10~ 3 cm~ 3 , which fixes Z to be 1.8 X 10~ 5 . Finally 
we take n c o(0)(= P c d(0)/ to p) = 180cm -3 and a c D = 2kpc (Malu- 
muth & Kirshner 1985) except for model E explained below for which 
ra cD(0) is varied. This corresponds to the velocity dispersion of 330 
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km/s from equation cr? D = 4:irGp cD (0)al D [3/2 - 2 In 2] (Tonry 1985), 
where <7 c rj is the velocity dispersion of stars in the cD galaxy. These 
parameters fixes the value of i? c D- Thus there is no free parameter in 
the source model(except for model E). 

Finally in the diffusion coefficient we take ay = lOOOkm/s where 
ay is the line-of-sight velocity dispersion or a r = a/y/3, and normalize 
%al(0) so that A500 = 21 as above. This makes _R e ff the control 
parameter of Do. 

For this parameter choice we are left with 5 free parameters, M c , 
(3, r s t, R e ff, and n c o(0). In order to see the effect of these parameters 
on the iron abundance distribution, we solve the evolution equation 
(18) for a series of models whose details are given in Table 1. 

Models labeled by A, B, C, D and E are ones to see the effect of 
cooling flows, for which D = 0. In the series A, B, C the first source 
model is used as S, while in the others the second source model is 
adopted, /e is put (e~' r / a + l)/2 in B models to see its effect but is 
put unity for the others. Models A and C are built to see the effect 
of (3 and r st . Since the fitted values of (3 are contained in the range 
0.5 < (3 < 0.7 for most clusters, we take 0.5 and 0.7 as typical values. 
Observationally p c B( r ) (r » « c d) of various clusters determined by 
velocity dispersions and core radii are in the range of factor 2 or 3 of 
the mean value in the 1 sigma level(Table 5 of Malumuth & Kirshner 
1985). Model E is to see the variation caused by this scattering. a c D 
is fixed to be 2kpc. 

On the other hand models F and G are taken to see the effect of 
galactic motions. Their model parameters are the same as those of 
model D except that the diffusion term is effective. The difference 
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between models F and G are the length of the period when the 
diffusive mixing is effective: = 10 10 yr = / i n model F and = 
f 8 yr = ^gw in model G. In both models we switch off the cooling flow 
to see the genuine effect of galactic motion. Further we set D(r) = 
for r < 17kpc for the following reason. Galaxies going through this 
region will be disrupted by tidal force of the cD galaxy. As a result 
these galaxies do not provoke mixing, but instead contribute to the 
source term, which is taken into account in the second source model. 
The tidal radius of a galaxy with mass m at distance r from the 
cluster center is expressed in terms of the cluster mass inside that 
radius, M(r) as 

R t = (m/SM(r)) 1/3 r (26) 

(Binney & Tremaine 1987). Since M(r) is dominated by the cD galaxy 
mass in the present case, for the typical values m = 5 X 1O 11 M and 
R t = fOkpc, this equation yields r ~ 17kpc. 

The value of Z at / = to in units of the solar abundance 4x 10~ 5 for 
these models obtained by numerical integration of Eq.(l) are shown in 
Figure 1 to 5. Further in Table 1 the two characteristic dimensionless 
velocities of the models, u a and m_r cD , and the abundance concentra- 



tion parameter Z defined as follows are shown: 

u(r) = v(r)t /r, (27) 

URcD = u(RcD), (28) 

u a = u(a), (29) 

Z = Z(0.1a,t )/Z(1.0a,t ). (30) 



Figure la, lb and Table 1 show that the cooling flow flattens the 
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initial iron inhomogeneity produced by the source term quite efficiently 
if the flow velocity is sufficiently large. The reason of this flattening is 
quite simple: it is an inward drift and compression of the distribution 
pattern by the cooling flow. To be precise, the value of Z is conserved 
in the Lagrange coordinate of the flow when the diffusion is negligible 
after the windy period. Since the iron abundance distribution just 
after the windy period becomes flatter in the outer region as r in- 
creases, the final distribution becomes flatter when the characteristic 
flow velocity u a increases, as is seen in Table 1. 

The cooling flow does not affect the value of Z at the center where 
M vanishes. This produces the sharp peaks at the center in Figures 
la, lb and 2a. For the same reason we get central peaks of width r st 
in model C as is clearly seen in Figure 2b. The inverse correlation 
between u a and Z for model C in Table 1 is due to these central 
peaks. If these peaks are observed in clusters with flat iron abundance 
distributions and strong cooling flows, they will strongly suggest that 
the flat distributions have been produced by cooling flows. However, 
as these peaks have narrow widths, they may not be observed by 
detectors with low spatial resolutions even if they exist. 

As is seen from comparison of Figures la, lb and 2a, the flattening 
effect of cooling flow becomes more prominent as the source term is 
more centrally concentrated for the source model with no cD contri- 
bution. In spite of this source term dependence Table 1 shows that 
the initial inhomogeneity produced by the source term is practically 
erased if u a is larger than unity. 

The abundance distribution just after the windy period in the 
source model with cD contribution is different from that in the source 
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model without cD contribution due to the existence of a broad, high 
central peak. This peak is made narrower by cooling flow as shown in 
Figure 3. As is seen in Figure 4, the final width of the peak is very 
sensitive to the value of n c o(0). This is because the initial width of the 
peak, i?cD 5 increases with ra c D(0) and u(r) decreases with r. Note that 
in spite of this sensitivity to physical parameters the phenomenologi- 
cal quantity w_r cD yields a good index to decide the fate of the peak. 
Roughly speaking, the peak disappears if the value of wr cD is around 
unity or greater. 

The final distributions in these models D and E resemble those in 
model C with r s t 7^ if the cooling flow is moderate. However, since 
the peak value is larger in the former than in the latter, they can be 
distinguished by observations with good spatial resolutions. 

Finally Figure 5 and Table 1 show that the efficiency of the diffusive 
mixing due to galaxy motion is sensitive to _R e ff an( i the duration of 
the diffusion. In particular from Figure 5a we see that the central peak 
due to the cD contribution is smeared out except in the region r < R t 
if i? e ff is larger than lOkpc during most of the cluster life. The cooling 
flow can easily erase this residual peak within R t if it reaches the 
cluster center. As is discussed in detail in the next section, however, 
it is likely that such strong diffusion does not occur in reality. Figure 
56 and Table 1 also show that the central peak by the cD contribution 
survives after the windy period even if _R e ff becomes as large as 30kpc 
due to galactic winds. 
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4 SPECIFIC MODELS OF THE OB- 
SERVED CLUSTERS 



As stated in the introduction, among the 12 clusters for which 
the iron abundance distribution has been observed, Virgo (Koyama et 
al. 1991), Centaurus (Fukazawa et al. 1994), AWM7 (Ohashi 1994; 
Fukazawa 1994), A496, A2142 and A2199 (White et al. 1994) have 
abundance gradients, but in Coma (Hughes et al. 1988, 1993), A1060, 
MKW3s (Ohashi 1994; Fukazawa 1994), Hydra-A (Ikebe et al. 1994), 
and A1795 (White et al. 1994), the gradients are small. Ponman 
et al (1990) and Kowalski et al (1993) reported the presence of the 
abundance gradient in the Perseus cluster. However, the recent ASCA 
observation gave a contrary result(Ohashi 1994). The physical param- 
eters of these clusters determined by observations are summarized in 
Table 2 where tip and are explained later. 

First from these data we construct specific models for each cluster 
within the framework given in §2. 

As the galaxy distribution function p ga ± we use the King model (6) 
for all clusters except for Virgo cluster and determine /O ga i(0) for each 
cluster with the help of Eq. (25). We assume the core radii of Eq. 
(6) are the same as those of Eq. (7) except for AWM7 and MKW3s. 
The cooling radii of AWM7 and MKW3s are determined from Figure 
1 of Canizares, Stewart & Fabian (1983). Since A500 is not known for 
MKW3s and Hydra-A, we tentatively take A 500 = 21 (Bahcall 1981) 
for these clusters. Finally for Virgo cluster, since the cooling flow 
should be regard as that of M87, we use the stellar mass density as 
^gal(r). 
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For most clusters p sas (r) is well fitted by the (3 model (7) outside 
the cooling radius, but for clusters with strong cooling flows the best 
fit (3 models yield smaller values in the central region than the observed 
values (Jones & Forman 1984). Since this deviation affects the cooling 
flow velocity, we modify the simple (3 model as follows for r < r c : 

PlJr) = (n d /n^- r ^p gas (r), (31) 

where tip is p g as(0)/ra p of the best fit (3 model p sas (r), and n& is the 
corresponding observed value. 

Further for AWM7 and MKW3s the gas density is not be fitted by 
the (3 model. Instead the following simple power-law model yields a 
good approximation(Kriss, Cioffi, & Canizares 1983): 

Pgas(r) = m p ra b (r/lkpc) Qf . (32) 

The best fit values of the parameters are n b = 0.32cm -3 and a = 
-1.04 for AWM7, and n h = 2.17cm" 3 and a = -1.35 for MKW3s. In 
these two clusters, the mass deposition rate is given at the radius of 
the central bin used in the observations, which are 70kpc for AWM7 
and 40kpc for MKW3s (Canizares, Stewart, & Fabian 1983). Taking 
account of these spatial resolutions of the observations, we modify 
Eq.(32) so that 

Pgas is constant within these central bins in order to 
avoid the divergence at the center. 

For all the clusters we assume to = 10 10 yr and ^gw = 10 8 yr. We 
also take f w = 0.18 and fp e = 6.8 X 10~ 4 as in the previous section. 

For Coma and A1060 we apply model A in the previous section 
because, though there are two cD galaxies in these clusters, they are 
not at rest at the centers. With f w and fp e given above, the density 
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distribution models explained above completely fixes the source term 
for these clusters. 

On the other hand for the other clusters we construct models by 
adding the diffusive mixing term by the galactic motion to model 
D. We use the fixed common values, a c D = 2kpc and /O c D(0)/m p = 
180cm -3 , to calculate the contribution from cD galaxies to the source 
term for all the clusters as we have no data on them. 

The study of Bahcall (1981) shows that there is no noticeable dif- 
ference in the line-of-sight velocity dispersion of galaxies between the 
clusters with centrally concentrated iron distributions and the others: 
its mean value is 854km/s for the former clusters, Virgo, Centaurus, 
A2142, A2199 and AWM7, and 960km/s for the latter, Coma, A1060, 
A1795 and Perseus. If the observed central concentrations are as- 
sumed to be produced by the central cD galaxies, this and the result 
for model F in the previous section suggest that i? e ff must be smaller 
than lOkpc around the center during most of the lives of the clusters, 
provided that the value of _R e ff does not vary from cluster to cluster. 
This is consistent with the facts that the elliptical fraction is high in 
the central region of clusters (e.g. Whitmore & Gilmore 1993) and 
elliptical galaxies have less gas than spiral galaxies, because a com- 
pletely stripped galaxy hardly disturbs intracluster gas as shown by 
Gaetz, Salpeter and Shaviv (1987) (see Fig. 2a in their paper). There- 
fore we take as _R e ff the value determined by the right hand of Eq.(12) 
for < / < ^gw an( i 5kpc for ^gw < t. This with the observed values 
of a and ra ga i determine the effective diffusion coefficient. As in the 
previous section ra ga i(0) is determined by A500 (Bahcall 1981) and we 
assume D(r) = for r < 17kpc. 
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The present profiles of the iron abundance distributions calculated 
for these specific models are shown in Figure 6(a) - 6(c). The corre- 
sponding characteristic velocities u a and m_r cD and the final density 
contrast index Z are listed in Table 3. Here Z is defined by Eq.(30) 
except for Virgo while the value defined by 

Z = Z(2',t )/Z(10 , ,t ) (33) 

is given for Virgo, assuming the distance to it is 15Mpc. The reason 
we take the values at 2' and 10' is that the iron abundance gradient is 
actually observed in this region by ASCA (Koyama private communi- 
cation). 

Figure 6(a) - 6(c) and Table 3 show that our model calculation 
gives iron abundance distributions which are consistent with the ob- 
servations for the clusters observed by ASCA (Coma, 1060, Virgo, 
Centaurus, AWM7, MKW3s, Perseus and Hydra-A) and A1795. This 
indicates that the combination of the strength of the cooling flows 
and the contribution of the central cD galaxies explains the variety 
of the iron abundance distributions for these clusters. In particular 
it can be understood that the cooling flows of Virgo, Centaurus and 
AWM7 have been too weak to modify the initial iron distributions 
while they have been erased by strong cooling flow compression in 
A1795, MKW3s, Perseus and Hydra-A (see u RcD of Table 3). 

On the other hand the situation for A496 and A2199 is subtle. 
Though our calculation gives central peaks for them, they are almost 
within 1'. Hence when they are averaged over 6' the angular scale 
corresponding to the field of view of Einstein SSS, they may not give 
as large central enhancement as observed. Finally in the case of A2142 
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the central peak disappears in our model contrary to the observation. 
As we mentioned in §3, however, there is an uncertainty in ra c rj(0) and 
the duration of the cooling flow, so the result is not conclusive. 

One interesting point of our results is that the uniform iron dis- 
tribution of Coma and A1060 are accounted for even if no merging 
has occured. Of course this does not imply that possible merging phe- 
nomena in the past are not relevant to the present iron abundance 
distribution. Instead it rather indicate that the contribution of the 
central cD galaxy is crucial in accounting for the present central con- 
centration of the iron in clusters with central cD galaxies because if 
the contribution were neglected, we would obtain flat distributions as 
in Coma and A1060. 

5 CONCLUSION 

In the present paper we have investigated the influence of cool- 
ing flow and galactic motion on the iron abundance distribution in 
clusters of galaxies by simple spherical models. We have performed 
numerical calculations for some typical models and shown that the 
cooling flow is very effective in flattening the centrally concentration 
of iron abundance produced during the initial windy period if the flow 
velocity is sufficiently large. We have also shown that the supply of the 
metal rich gas by the central cD galaxies are essential in accounting 
for the existence of clusters with highly central concentrated distri- 
bution of iron abundance, and that the diffusive mixing by galactic 
motion should be small enough for those central peaks to survive. The 
latter is consistent with the fact that the central regions of clusters 
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are usually dominated by gas-poor elliptical galaxies, but gives the 
constraint that the effective radius of galaxies should be smaller than 
30kpc during the windy period. 

Further on the basis of this general analysis we have constructed 
specific models for each of the observed clusters. We have shown that 
the observed variety in the iron abundance distribution of clusters with 
central cD galaxies at rest at center can be qualitatively explained in 
our scheme as a result of the production of a centrally concentrated 
iron distribution by the cD galaxy and its subsequent deformation by 
cooling flow, except for a few clusters for which the available data 
have too poor resolution or some important physical parameter are 
unknown. This result is remarkable taking into account the simplic- 
ity of our model, and we believe that it suggests the validity of our 
scenario. 

One important feature of our model is the existence of a clear ob- 
servable prediction. Since the flat distribution of the iron abundance is 
produced by the central compression of the iron rich region toward the 
center in our model, a sharp peak in the iron abundance distribution is 
left at the center if the cooling flow does not reach the center. If such 
a peak is detected by future observation with high spatial resolutions, 
it will give a strong evidence for the validity of our scenario. 

Of course such a peak may not be left if the effective radii of 
galaxies become larger than 30kpc during the windy period. If we 
extend our argument to other heavy elements with this constraint 
in mind, we obtain another interesting prediction. According to the 
study of the chemical evolution of galaxies (e.g. David et al. 1991b) 
strong galactic winds are produced by Type II supernovae which eject 
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oxygen rich gas, while iron enriched gas is ejected by Type I supernovae 
at low ejection rate after the strong wind period. Thus, if our scenario 
is correct, this leads to the prediction that the iron distribution has 
a central peak but the oxygen abundance is uniform in clusters of 
galaxies which have central cD galaxies at the centers and have weak 
cooling flows. The observational check of the prediction will give useful 
information on the structure of galaxies during the windy period as 
well as their chemical evolution. 

We would like to thank K.Koyama, A.Habe and T.Tsuru for useful 
suggestions and providing us new data of ASCA. We also wish to 
thank K.Tomita, J.Yokoyama, Y.Yamada, R.Nishi, Y.Fujiwara and 
N.Makino for discussion. 
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Figure Captions 



Fig.l.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center, (a) model Al. (b) model A2. 
Fig. 2.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center, (a) model B. (b) model C. 
Fig. 3.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center for model D. 
Fig. 4.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center for model E. 
Fig. 5.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center, (a) model F. (b) model G. 
Fig. 6.- Iron abundance distributions in units of the solar abundance as a function 

of the distance from the cluster center, (a) Coma, A1060, Virgo and Cen- 

taurus. (b) A496, A2142, A2199 and AWM7. (c) A1795, MKW3s, Perseus 

and Hydra-A. 
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TABLE 2 
OBSERVATIONAL DATA 
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Note: In AWM7 and MKW3s M c is the mass deposition rate at the innermost 
observational bin. In A496 and A2142, (3 is assumed to be 1. We used 
H = 50km s" 1 Mpc" 1 . 

a. Edge k Stewart 1991a, b. Takano 1990, c. Matilski et al. 1985, d. Kriss 



et al. 1983, e. Canizares et al. 1983, /. Hughes et al. 1988, g. Jones 
k Forman 1984, h. David et al. 1990, i. Abramopoulos k Ku 1983, j. 
Fabricant k Gorenstein 1983, k. Zabludoff et al. 1990, /. Binggeli et 
al. 1987, to. Lucey et al. 1986, n. Hintzen k Scot 1979, o. Gregory k 
Thompson 1984, p. Materne k Hopp 1983. 
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urcD, u a AND Z OF EACH CLUSTER 



Cluster 


URcD 




Z 


Coma 




0.001 


1.3 


A1060 




0.1 


1.1 


Virgo 




1.5 


(3.9) 


Centaurus 


0.6 


0.3 


4.4 


A496 


0.6 


0.2 


3.8 


A2142 


1.0 


0.03 


1.1 


A2199 


0.7 


0.3 


2.2 


AWM7 


0.5 


0.1 


2.7 


A1795 


1.4 


0.2 


1.2 


MKW3s 


2.5 


0.3 


1.0 


Perseus 


1.3 


0.2 


1.0 


Hydra-A 


3.2 


1.7 


1.3 



